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Abstract 

We study dynamics of a probe particle, which performs biased diffusive motion in a one- 
dimensional adsorbed monolayer of mobile hard-core particles undergoing continuous exchanges 
with a vapor phase. In terms of a mean-field-type approach, based on the decoupling of the third- 
order correlation functions into a product of the pairwise correlations, we determine analytically 
the density profiles of the monolayer particles, as seen from the stationary moving probe, and 
calculate the terminal velocity V pr , mobility fi pr and the self-diffusion coefficient D pr of the probe. 
Our analytical results are confirmed by Monte Carlo simulations. 



Key words: Hard-core lattice gas, Langmuir adsorption/desorption model, tracer diffusion and 

mobility. 



1 Introduction 



When a solid surface is brought in contact with an ambient gas phase, the interactions between the 
molecules of both systems often result in formation of a layer of gas particles covering the solid surface. 
Following a seminal work of Langmuir (see, e.g., in Ref.Jl|), who has invented a somewhat simplified 
description in which the interactions between the adsorbed molecules were regarded as a mere hard- 
core, thermodynamic properties of such adsorbed layers, as well as different forms of possible phase 
transformations have been extensively studied and a number of important developments have been 
made. In particular, subsequent studies included more realistic forms of intermolecular interactions 
or allowed for the possibility of multilayer formation. In consequence, more sophisticated forms of 
adsorption isotherms have been evaluated which explain quite well available experimental data (see, 
e.g. Refs. 11,0). 

Considerable effort has been also invested in understanding of molecular diffusion in adsorbed lay- 
ers, which has a strong impact on their properties @, ||. Here, some approximate results have been 
obtained for both dynamics of an isolated adatom on a corrugated surface and collective diffusion, de- 
scribing spreading of the macroscopic density fluctuations in interacting adsorbates being in contact 
with the vapor phase [f§, H, H 0, H II- O n the other hand, available studies of the tracer diffusion 
in adsorbed layers, which is observed experimentally in STM or field ion measurements and provides 
a useful information about such properties of monolayers as, e.g., their intrinsic viscosity, pertain to 
stricktly two-dimensional models excluding the possibility of particles adsoption or desorption (see, e.g. 
Refs.[||], [|, and references therein). Analysis of the tracer diffusion and mobility of impure molecules 
in adsorbed monolayers undergoing exchanges with the vapor phase seems to be lacking at present. 

In this paper we study the time evolution of a model system consisting of a solid substrate covered 
by a monolayer of mobile hard-core particles undergoing continuous exchanges with a vapor, and a 
single impure, probe particle, which is subject to a constant external force E and hence performs 
a biased random walk constrained by hard-core interactions with the monolayer particles. Here, we 
concentrate on the one- dimensional case and model, in a usual fashion, the solid substrate as a regular 
one-dimensional lattice of adsorbing sites, which can support, at most, a single occupancy; results for 
the two-dimensional monolayer will be published elsewhere. In terms of a mean-field-type approach 



of Ref.jll], which pressumes certain decoupling of the third-order correlation functions, we define the 
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density profiles of the monolayer particles, as seen from the stationary moving probe, and determine 
analytically the probe terminal velocity V pr (E). In the most general case V pr (E) is obtained implicitly, as 
a solution of a non-linear equation relating its value to the system parameters. This equation simplifies 
in the limit of small E and here the probe velocity can be found explicitly; we show that it obeys 
Vpr(E) ~ E/(, where the friction coefficient (, which is inverse of the probe mobility fi pr , is expressed 
through the microscopic parameters characterizing the system under study. This result establishes the 
frictional drag force exerted on the probe by the monolayer particles in the \ow-E limit and thus can be 
thought off as the analog of the Stokes' law for the one-dimensional monolayer in contact with vapor. 
Lastly, we determine the self-diffusion coefficient D pr of the probe, which is computed here by assuming 
heuristically the validity of the Einstein relation between the self-diffusion coefficient and the mobility. 
Our analytical results for the probe terminal velocity and the self-diffusion coefficient, as well as for 
the stationary density profiles around it are confirmed by Monte Carlo simulations of the corresponding 
master equation by the method of Gillespie ||15|| . 

The paper is structured as follows: In Section 2 we formulate the model and introduce basic notations. 
In Section 3 we write down the dynamical equations which govern the time evolution of the monolayer 
particles and of the probe. Sections 4 and 5 are devoted to the analytical solutions of these evolution 
equations in the continuous-space limit and on the discrete lattice, respectively. In these sections we also 
present a comparison of the analytical results and Monte Carlo simulations data. Finally, we conclude 
in Section 6 with a brief summary and discussion of our results. 

2 The model 

The model consists of a one- dimensional regular, infinite in both directions lattice of spacing a, the 
sites X of which are connected to a reservoir (a vapor phase) containing an infinite number of identical, 
electrically neutral gas particles (Fig.l). The particles from the reservoir may adsorb onto the lattice 
sites, desorb from them or move along the lattice by performing symmetric random walk between the 
neighboring sites; adsorption onto the lattice and hops between the lattice sites are constrained by 
hard-core exclusion - each lattice site can be either singly occupied or vacant. The state of each lattice 
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site X at time t is described by time- dependent occupation variable n(X), which can take two values 



As one may readily notice, the just- described system corresponds to a one- dimensional version of the 
Langmuir adsorption/desorption model, with the only difference being that the lateral diffusion of the 
adsorbed molecules is allowed. Note, however, that the possibility of particles lateral diffusion makes 
the model more complicated, compared the Langmuir's one, since now the evolution of particle local 
densities at different sites is coupled due to diffusion. 

Further on, we place at the site X = at time t = an extra particle, which will be referred to in 
what follows as the "probe", since it allows us to probe the resistance offered by the adsorbed layer to 
an external perturbance. Position of this particle at time t for a given realization of the process will be 
denoted as X pr (t). We stipulate that this particle is different from the monolayer particles in that it 
can not desorb back into the vapor phase, i.e. is constrained to move along the lattice only. Second, we 
suppose that this only particle is charged (for simplicity, we set the charge equal to unity in what follows) 
and is subject to a constant external electric field E, which favors its motion in a preferential direction. 
The questions which we address here are, first, the dependence of the probe particle terminal velocity 
Vp r (E) = limt-tcaVprit) on the magnitude of the driving force E and other system parameters; second, 
the form of the density profiles as seen from the stationary moving probe, and lastly, the self-diffusion 
coefficient of the probe in absence of the driving force. 

Now, we define particle dynamics and system parameters more precisely: 

(a) Diffusion and desorption. Each of the adsorbed at time moment t gas particles waits a random, 
exponentially distributed time with mean r*, and then selects between either of three possibilities; it 
may choose to leave the lattice with a probability g, or attempt to hop, with the probability (f — g)/2, 
to one of the two neighboring sites. If the desorption event is chosen - the particle leaves the lattice 
instantaneously. On the other hand, in case when the particle attempts to hop to one of the neighboring 
sites, the jump can actually occur only if the target site is empty; otherwise, the particle remains at its 
position. 

(b) Adsorption. Particles of the reservoir wait a random, exponentially distributed time with meanF] 
1 We choose for simplicity the same mean waiting time for the particles of the lattice and those in the reservoir; our 



1, if the site X is occupied by a gas particle 




0, otherwise 
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T a d = t* and then attempt to adsorb onto the lattice with a probability f. As previously, an adsorption 
event takes place only if the chosen site is empty 

Note, that the total number of particles on the lattice is not conserved in such a dynamics. Mean 
density p =< r)(X) >, however, approaches as t — > oo a constant value 

" = 7T? (2) 

which relation is often called in the literature as the Langmuir adsorption isotherm |1] . Here we suppose 
that parameters / and g are independent of each other and may take arbitrary values from the interval 
[0; 1]. In reality, their values are prescribed by the pressure of the gas phase, precise form of the solid- 
gas interaction potential and the temperature. As well, the latter determine the characteristic time t*, 
which is also considered here as an independent given parameter. 

(c) Dynamics of the probe. The probe particle waits an exponentially distributed time with mean r, 
(which can be, in general case, different of r*) and then selects, at random, a jump direction: It chooses 
a right-hand or left-hand adjacent site with probabilities p and q = 1 — p, respectively. Similarly to the 
gas particles, the jump is only then fulfilled when the selected site is vacant at this moment of time. 

In a usual fashion, the jump probabilities are related to the external electric field E, (oriented in the 
positive direction, such that E > 0), and the reciprocal temperature (3 by 

p/q = exp(flo-E), (3) 

where we have set, for simplicity, the probe charge equal to unity. This relation together with the 
condition p + q = 1 defines the values of p and q. 



3 Evolution equations 

Let P(X pr , rj; t) denote the probability of finding at time t the probe at the site X pr and all other 
adsorbed particles in the configuration 77 = {i](X)}. Further on, let r/ 2,2+1 be the configuration obtained 
from 77 by exchanging the occupation variables of sites X = z and X = z + 1, i.e., r](z) <-> r\{z + 1), 
and if - the configuration obtained from r\ by replacement r\{X = z) — > 1 — r](X = z). Then, summing 
results can be simply generalized for arbitrary T a d by replacement / — > r* / /r ac i. 
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up all events which may change or result in a configuration (X pr , 77), we find that the time evolution of 
P(X pr , 77; t) is guided by the following master equation: 

P(X pr , m t) = ^ £ {P(X^>* + V)-P(*pr, 

Z^Xpr — (J ^Xpr 

+ ^{(1 - V (X pr ))P(X pr - a, 77; t) - (1 - V (X pr + a))P{X pr , 77; £)} 
+ i {(1 - 77(X pr ))P(X pr + a, 77; f) - (1 - 77(X pr - a))P(X pr , 77; /;)} 
+ ^ E {(l-^))P(X pr ,77 2 ;t)-^)P(X pr ,77;t)} 

2^ Xpr 

+ ^ E {^)P(X pr ,^;t)-(l-7 / (z))P(X pr ,7 ? ;t)} (4) 

where the dot denotes the time derivative, the terms in the first three lines describe respectively the 
diffusion of the adsorbed gas particles and of the probe, which proceed due to the Kawasaki-type particle- 
vacancy exchanges, while the last two lines account respectively for the Glauber-type desorption and 
adsorption processes. 

Equation ([|) allows to compute the instantaneous velocity of the probe molecule. Multiplying both 
sides of Eq.(B) by X pr , and summing over all possible states (X pr , 77) of the system, we find that : 

V pr (t) = = a -{p{l - k(a; t)) - q(l - k(-a; t))} (5) 

where X pr (?;) denotes the mean displacement of the probe at time t and 

k(X;t) = J2 rj{X pr + \)P(X pr ,rj;t) (6) 

(Xp r , ri) 

defines the probability of having at time t an adsorbed gas particle at distance A from the probe, or in 
other words, can be interpreted as the monolayer density profile as seen from the moving probe. 

Hence, in order to compute V pr (t) we have to determine the evolution of k(X;t). From Eq.(£|) it 
follows that we have to consider separately the evolution of k(X;t) for |A| > o and |A| = a, since the 
monolayer particles dynamics is different in these two domains. We start first with the case |A| > a, in 
which domain the evolution of k(X; t) is not directly affected by the presence of the probe. In this case 
we find from Eq.(^) the following equation 

HA; t) = ^^{k(\ + a;t) + k(\-a;t)-2k(\;t)}-t±Ak(\;t) + ^- + 

+ - E U ~ v(X pr + cr))P(X pr , n; t)^7](X pr + A + cr) - 77(X pr + A)} - 

(X pr ,v) 

~ \ E (l-v(X pr -cr))P(X pr , m t){ V (X pr + X)-r ] (X pr + X-a)}, (7) 

(X pr ,v) 
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where the terms in the first line represent the contributions due to the particles diffusion, desorption and 
adsorption, while the terms in the last two lines are associated with the translation of the configuration 
r) — > 7] z ^ z±l due to the backward and forward hops of the probe particle. These terms are non-linear with 
respect to the occupation variables and thus couple the evolution of the "pairwise" correlation function 
k(X; t) to the evolution of the third-order correlations. Thus, to define the behavior of k(X; t) one 
faces the problem of solving an infinite hierarchy of coupled differential equations for the higher-order 
correlation functions. 

Here we will resort to an approximate, mean-field-type decoupling scheme, which has been first 
applied in Ref. to describe tracer diffusion in a one- dimensional hard-core lattice gas with conserved 



number of particles, i.e. a gas for which the exchanges with the reservoir are forbidden and both 
/ and g are equal to zero. It has been shown in Ref.[|ll[] that results of such an approach are in a 
very good agreement with the Monte Carlo simulations data. Moreover, rigorous probabilistic analysis 
HJ of this model produced essentially the same results as Ref. pif , thus proving that such a decoupling 



scheme renders exact description of the model. We note also parenthetically that a very good agreement 
between the numerical data and analytical predictions, based on the same decoupling scheme, have been 
observed for a slightly different model of hard-core gas spreading on a one- dimensional lattice from a 



reservoir connected to one of the lattice sites |jl2 |. We thus adopt here this mean- field-type approach, 
which provides quite a good description of related dynamical models, and will verify in what follows 
our analytical predictions against the results of Monte Carlo simulations. 



The decoupling scheme of Ref. jTIJ is based on the assumption that the average with the weight 
P(X pr , 77; t) of the product of several occupation variables of different sites factorizes into the product 
of their average values with the weight P(X pr , r\\ t). Namely, it assumes that the average product of the 
occupation variables factorizes as 

J2 v(X pr + A)(l - V (X pr ± a))P(X pr , 77; t) = 
= { ]T r]{X pr + X)P{X pr ,r];t)}x{ £ (1 - r){X pr ± a))P{X pr , m t)} = 

(X pr ,ri) (X pr ,ri) 

= k(X;t)(l-k(±a;t)) (8) 
Then, taking advantage of Eq.(|8|), we can rewrite Eq.(^), which holds for |A| > a, in the following form: 

k(X;t) = ] -^{k{X + a-t) + k(X-(r-t)-2k{X-t)}-t^-k{X-t) + l- + 

2r* t* t* 
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+ £{1 - k(a; t)}{k(X + <T]t)- fc(A; t)} - -{1 - k(-a; t)}{k(X; t) - k(X - a; t)} (9) 
r r 

Note, however, that despite the fact that the decoupling in Eq.(|3|) allows us to close the hierarchy in 
Eq.(0) at the level of pairwise correlations, the resulting equations still pose some technical problems 
for solving them; namely, they are non-linear with respect to k(X; t) differential equations. The method 
of solution will be discussed in the next two sections. 

Now, similar analysis can be carried out to derive the dynamical equations, which govern the evolu- 
tion of k(X;t) at points |A| = a. We find, respectively, 

fat) = ^-{k(2a;t)-k(a;t)}-^-k(a;t) + ^ + 

+ -{-qk(a;t)(l-k(-a;t))+p(l-k(a;t))k(2a;t)} (10) 



and 



k(-a;t) = l -^L{k{-2o-t) - k(-a;t)} - Ltl k (- a . t ) + L 



+ -{-pk(-a;t)(l-k(a;t)) + q(l-k(-a;t))k(-2a;t)} (11) 

Equations ©,© and ([lO| ) , ([lT|) constitute a closed system of non-linear equations, which suffice the 
computation of the probe velocity and other characteristic properties. 



4 Stationary solution of the evolution equations in the conti- 
nuous-space limit. 

Consider first the solution to Eqs.©,@ and (|10|), (|ll]) in the continuous-space limit. Expanding k(X± 
a;t) in Taylor series up to the second order in powers of o (diffusion limit), we have that k(X;t) with 
|A| > obeys: 



H X; t)^D ^> + V At) d -^> - l±l H X; t) + l, (!2) 

where D denotes the "bare" diffusion coefficient of the adsorbed gas particles, D = (1 — g)a 2 /2t* , and 
the velocity V pr (t) is now given by 



V pr (t) = ^-{l - k(X = +0; £)} - ^{l - k(X = -0; t)} (13) 
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Similarly, we find that Eqs.(|9|) and (|H]) take the form 

dk(X;t) 



and 



ak(X = +0;t) = D 



ak(X = -0;t) = - A 



dk(X;t) 



A=+0 



+ {v pr (t) - (/ + 9)^}k(X = +0; t) + ^ 



°1 

T* 



dX 



A=-0 



~{V pr (t) + (f + g)^}k(X = -0;t) + 



of 



(14) 



(15) 



We turn next to the limit t — > oo. Assuming first that the probe terminal velocity approaches 
some constant value V pr (E) as t — > oo, i.e. V pr (E) = /im t ^ 00 V^ r (t), we find the stationary solution 
fe(A) = limt^k^t) of Eqs.(0),(0) and ([15]). This reads: 

(16) 



k(X) = p s 1 + A± exp ( - |A|/A± 



In Eq . (|i~6"D p s is determined by Eq.(|), the sign "+' 
lengths A± are given by 



') corresponds to A > (A < 0), the characteristic 



A- 



\ 



V pr (E) a^2 ± Vpr(E) (T-,-1 



2A, 



2Dn 



(17) 



where the parameter cr e // has the dimensionality of length, cr e // = D r* /a(f + g), while the amplitudes 
At in Eq.(||) obey 



A- 



± 



(18) 



D n + a e// /A T 

Now several comments on the results in Eqs.fll~6| ) , (lT7l) and (]T8|) are in order. Note, first, that A_ > A + , 
and consequently, the local density past the probe approaches its non-perturbed value p s slower than 
in front of it; this signifies that correlations between the probe position and particle distribution are 
stronger past the probe and demonstrates the memory effects of the medium. Next, A + is always 
positive, while v4_ < 0; this means that the density profile is a non-monotoneous function of A and is 
characterized by a jammed region in front of the probe, in which the local density is higher than p s , and 
a depleted region past the probe in which the density is lower than p s . As a matter of fact, the depleted 
region appears to be more pronounced than the jammed one and the mean density in the adsorbed 
monolayer perturbed by the probe is lower than the mean density in the unperturbed monolayer; one 
readily finds that the integral deviation Q of the density from the equilibrium value p s , i.e., 

^ dX 



r — (*w 



Pi 



(19) 
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is always negative, 

n = -p a ( v ^ E)a '") 2 + l^ + (VAE)* e 7Ky (20) 

which means that the perturbance of the monolayer created by the probe shifts the balance between 
adsorption and desorption towards particles desorption into the reservoir. This result may be a bit 
counterintuitive and one could expect the reverse sign of Q, since the local density in front of the probe 
is higher than p s and particle depletion past the probe favors adsorption from the reservoir. However, 
Eq . (|20|) states that the reservoir fails to equilibrate the depleted region past the probe, whose size A_ 
exceeds the size A + of the jammed region. Apparently, this feature may be specific to a one- dimensional 
system only and the sign of Q may be different for 2D, in which case k(X) shows non-monotoneous 



A-dependence for A < in case of a lattice gas with conserved particles number JT0|- Note, lastly, that 
Q in Eq. (p0|) is proportional to the second power of the probe terminal velocity at small Peclet-type 
numbers P = V pr (E)a e ff/D , and increases linearly with V pr (E) when P is large. 

Now, we are in position to calculate the terminal velocity of the probe molecule. Substituting Eq . (|16|) 
into Eq. (|13D , we arrive at a closed-formed non-linear equation of the form 



(p - q)(l - p s )a V pr (E)a ef f 



P 



= r ~ P'° nr 1+/ /A + 1+„ l\ (21) 

r D q t L i + o- e ///A_ l + <T e ///A + J 

which determines V pr (E) implicitly. In Eq.([H]) the first term on the rhs is a trivial, mean-field-type 
result, which obtains in the limit of, say, perfectly stirred monolayer with D — > oo (or r* — > 0). 
The second term on the rhs of Eq . (]2l|) , which is proportional to the Peclet-type number P, has a more 
complicated origin and is associated with collective effects - formation of a non-homogeneous, stationary 
density profile around the probe, whose characteristics depend themselves on the velocity. 

Eq.(^l|) is rather complicated and can be solved for arbitrary values of the system parameters only 
numerically. It simplifies considerably, however, in the limit E — > 0, in which case V pr (E) can be 
calculated analytically. Expanding p, q and V pr (E) in powers of E and retaining only linear in E terms, 
we find that in this limit V pr (E) attains the following, physically revealing form 

V pr (E) « C^, (22) 



which can be thought off as the analog of the Stokes formula for driven motion in a one-dimensional 
monolayer undergoing continuous exchanges with the vapor phase. The friction coefficient in Eq.(p2|) is 

9 



given explicitly by 

C - 2r [n a P^ff 1 (23) 

(3a 2 (l-p s )L D r(l + ^a eff /ay 

Note, that ( is a sum of two contributions. The first one, ( m f, has an essentially mean-field-type form 
and is just a usual expression for the inverse mobility of a particle performing biased random walk on a 
one-dimensional lattice, divided by the fraction of non-occupied lattice sites, (1 — p s ), which is the mean 
density of vacant sites in a monolayer with homogeneous particle distribution, i.e. monolayer being in 
equilibrium with the vapor phase. On contrary, the second term, ( coop , stems out of the cooperative 
effects, associated with formation of a stationary, non-homogeneous density profile around the probe, 
and represents the net resistance offered by the monolayer particles to the probe. It may be instructive 
to single out this contribution explicitly; setting r = 0, (which means physically that the probe slides 
on the surface without dissipation regardless of the surface corrugation), we have 

C = Kil \ (24) 

coop (3a\l - p s ) (/ + g){1 + ^D r*/a 2 (f + g)) 

Comparing ( m f and Ccoopi we infer that the latter becomes progressively more important in the limit, 

for instance, when both / and g tend to zero (while their ratio f/g is kept constant, f/g — const = 

(1 — p s ) / p s ). Of course, such a behavior can be expected on physical grounds since this limit corresponds 

to a lattice gas with a fixed number of particles and suppressed exchanges with the vapor phase; in this 

case the friction coefficient (as well as the characteristic lengths A±) is no longer constant, but rather 



diverges as ( ~ X pr (t) ~ t 1 ^ 2 in the limit t — > oo [11, 13 



Consider finally the situation with E = 0, in which case the terminal velocity vanishes and one 
expects conventional diffusive motion with the mean square displacement of the form 



X*.(t) = 2D pr t, (25) 

where D pr is some unknown function of the system parameters. Heuristically, we can compute D pr 
for the system under study if we assume the validity of the Einstein relation p pr = f3D pr between the 
mobility p pr , p pr = UmE->o(V pr (E) / E) = and the self- diffusion coefficient D pr of the probe particle 



(see for more details 10). We find then that the probe self-diffusion coefficient should be given by 



_ a\l-p. 
pr ~ 2r 



(26) 



D t(1 + ^Jcr eff /a)- 

In order to check our analytical predictions in Egs-flOp, (f2l|) and fl26[) , we have performed numerical 
Monte Carlo simulations of the system evolution using the method of Gillespie [0. Results of these 
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simulations, performed at different values of the parameters /, g and p, are presented in Figs. 2 to 
4. These figures show that there is some systematic error between the analytical predictions of this 



Section, Eqs.(|T6|), (2j|) and (|26|) , and the numerical data; this error seems to be small when / is small 
but increases with increase of / or decrease of g, which means, apparently, that the error grows with 
the mean density of the monolayer. In particular, the maximal relative error R appears for g = 0.3 and 
amounts to 5.3, 7.7 and 7.9 per cent for Figs. 3,4 and 5, respectively. Hence, the agreement is quite fair, 
and consequently, Eqs.(|T6|),(|21"|) and ( p6| ) can be regarded as rather accurate approximate results. 

Now, one may attribute the origin of the discrepancy between the numerical and analytical results 
to the following two reasons. The first is, evidently, the decoupling of the hierarchy of equations for 
the correlation functions in Eq.(|S|), which allowed us to close this hierarchy at the level of pairwise 
correlations. Despite the fact that such a decoupling scheme yields exact results for tracer diffusion in 
one-dimensional hard-core lattice gases with conserved number of particles Jll], |13 |, it may well be that 
it fails for systems in which the total number of particles is not conserved. Second, such a discrepancy 
may result from the fact that here solution of the discrete-space equations has been found by turning 
to a continuous-space limit. As we proceed to show in the next Section, the latter point appears to be 
the most decisive and the solution of the discrete-space equations (|5|) , (p[) , fllCf) and ([LI]) is in a very good 
agreement with the numerical data. This implies, in turn, that the decoupling in Eq.(^) is also quite a 
plausible assumption for the system under study. 



5 Stationary solution of the discrete-space evolution equa- 
tions. 

Consider now the solution of the discrete-space (§) , (|10D and (|IID in the limit t — > oo. Denoting k n = 
k(X = na), and introducing auxiliary variables 

Zl = l + ^(l-h) (27) 

and 

= l + (28) 
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in order to avoid too lengthy formulae, we can rewrite Eq. (||) in form of the following recursion relation 

(29) 



ZiK+i - [zi + z_i + ]k n + z-xk n -i + p s -^— = 0, 

a eff °eff 



which holds for \n\ > 1. Equation (|29"D has to be solved subject to the boundary conditions 

z x k 2 - [z-i + — \h + p s — = 0, (30) 

°"e// CTeff 

and 



2;_ifc_ 2 - H + p s 







(31) 



Solution to these equations can be conveniently obtained by applying the discrete-space Fourier trans- 
formation, which yields an equation of essentially the same form as Eq. (|T6|) , i.e., 

k n = Ps 1 + A' ± exp{ — a\n\/X' ± 



(32) 

with, however, different amplitudes A'± and different characteristic lengths A' ± compared to those given 
by Eqs.(|18D and ([17]); we use here the prime to distinguish between the discrete and the continuous-space 
solutions. In the discrete-space Eq.(^) the characteristic lengths obey 



-i 



A' ± = =F o In 
while the amplitudes are given respectively by 

and 

A' 



zi + Z-i + (er/er e //) =F \ Ui + Z-i + (a /a, 



J eff, 



4zi^_i 



2zx 



(33) 



z x - Z_i 



Z-i — Ziexp(—a/X' + ) 



z\ - Z-i 



(34) 



(35) 



z_iexp(— er/A'_) — Z\ 

Now, we are in position to obtain a system of equations, determining implicitly the unknown pa- 
rameters Z\ and which will allow us to compute the terminal velocity of the probe molecule; the 
latter is related to z±\ by V^.{E) = D (zi — z_i)/<7. Substituting Eq.(|32"D into Eqs.(f27|) and (f28|), we 
find 



z\ = 1 + 



per 



D t 



1- ps- p s 



Z\ - Z-x 



Z-\exp(a / X' + ) — Zi 



and 



z-i = 1 + 



c/cr 



1 - Ps - Ps 



12 



z\ - 2-1 



2_i — z\exp(a/\'_ 



(36) 



(37) 



For arbitrary values of p, f and g the parameters z±i, defined by Eqs.(|36|) and fl3"T|), and consequently, 
the terminal velocity V'(E) can be determined only numerically (see Figs. 3 and 4). However, Vp r (E) 
can be found analytically in the explicit form in the limit of a vanishingly small force E, E — > 0. 
Expanding z±\ in the Taylor series in powers of E and retaining only linear with E terms, i.e. setting 



z±i ~ a + aj±iE, 



(3J 



where a = 1 
equations: 

and 



<t 2 (1 — p s )/2D r, we find that coefficients 7±i are determined by the system of two linear 



"7-i 



cr 2 p s 



(71 - 7_i) 



2D r 



exp[a/X' + (E = 0) 
(7i - 7-i) 



+ 



/3a 3 (l-p s ) 

ADnT 



(3cr\l-p s ) 



2D r 



exp 



(a/\'_{E = 0)) 



4ZV 



which yield the Stokes- type law in Eq.( 

2r 



with the friction coefficient 
apsVeff 2 



P<T*(1-Ps 



1 + 



(39) 



(40) 



(41) 



D or 1 + yJl + Aaaeff/cr 1 
Note, that similarly to Eq. (p3|) the friction coefficient calculated from the discrete-space evolution 
equations, is a sum two terms; the first one describes a trivial, mean-field-type behavior and coincides 
with the result obtained within the continuous-space approximation. The second one, which is associated 
with the formation of a stationary, non-homogeneous density profile around the probe, has a different 
form compared to that in Eq. ( |23| ) and reduces to it only in the limit a a e f f/a 1 . 

Finally, Eq. (|4*T|) allows to obtain the probe self-diffusion coefficient explicitly. Assuming again the 
validity of the Einstein relation, we find 

2 (l - P„) r. VPs&eff 2 1-1 



Dl 



pr 



2r 



1 



(42) 



d ot 1 + yjl +4aa eff /or- 
which has a different form compared to Eq. (p6|) . 

Comparison between the discrete-space solution and numerical data is presented in Figs. 2 to 5. We 



see now that Eqs.(|36),(|37),(|32) and (|42 ) agree with the Monte-Carlo results essentially better than their 
continuous-space counterparts; the maximal relative error is now independent of g and /, and amounts 
to only 1.9, 0.8 and 3 per cent for Figs. 3,4 and 5, respectively. We note finally that despite this good 
agreement we certainly can not claim that Eqs.(|36|),(^7|),(^) and (f42|) provide an exact solution of the 
model, which can be approached apparently by the method developed in Ref.|fL3"|. 
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6 Conclusions. 



To conclude, we have studied dynamics of a driven probe molecule in a one-dimensional adsorbed 
monolayer composed of mobile, hard-core particles undergoing continuous exchanges with the vapor 



phase. Within the framework of a decoupling procedure of Ref . [|TT| , based on the decomposition of 
the third-order correlation functions into a product of pairwise correlations, we have derived dynamical 
discrete-space equations describing evolution of the density profiles, as seen from the moving probe, 
and its velocity V pr (E). These equations have been solved both in the continuous-space (diffusion) 
limit, as well as in their original discrete-space form, which allowed us to check the appropriance of the 
continuous-space description. In both cases we have determined the probe particle terminal velocity 
V pr (E) implicitly, as a solution of non-linear equations relating its value to the system parameters. 
We have shown that the discrete-space solution provides a very good agreement between analytical 
and numerical results, while the continuous-space approach can be regarded only as a rather accurate 
approximation. Further on, we have found that in the limit of a vanishingly small driving force the 
probe velocity attains the form V pr (E) « -B/C, where E is the driving force and the friction coefficient 
( is expressed through the microscopic parameters characterizing the system under study. This result 
establishes the frictional drag force exerted on the probe by the monolayer particles in the low-E limit 
and thus can be thought off as the analog of the Stokes' law for the one-dimensional monolayer in 
contact with the vapor phase. Lastly, we have determined explicitly the self- diffusion coefficient D pr of 
the probe, which result has been also confirmed by Monte Carlo simulations. 



Acknowledgments. 

The authors wish to thank G.Tarjus, J.Talbot and J.Lyklema for helpful discussions. This work was 
supported in part by the French-German collaborative research program PROCOPE. 



7 Figure Captions 

Fig.l. One-dimensional lattice partially occupied by identical hard-core particles (filled circles) under- 
going exchanges with the reservoir - the vapor phase, g, f and I, I = (1 — g)/2, denote respectively 
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particles desorption, adsorption and hopping probabilities. The open circle denotes the probe molecule, 
whose hopping probabilities are p and q, respectively. 

Fig. 2. Density profile around stationary moving probe molecule for / = 0.1, g — 0.3 and p = 0.98. 
The solid line is the plot of the solution, Eq.(^), of the discrete-space evolution equations. The empty 
triangles denote the corresponding solution in the continuous-space limit. Filled squares are the results 
of Monte-Carlo simulations. 

Fig. 3. Terminal velocity of the probe molecule as a function of the adsorption probability / at 
different values of the parameter g. The probe hopping probabilities are p = 0.6 and q = 0.4. The 
solid lines give the solution of Eqs.(^) and (0), the dashed lines - of Eq. (|21~D , while the filled squares 
denote the results of Monte-Carlo simulations. Upper curves correspond to g — 0.8, the intermediate - 
to g = 0.5 and the lower - to g = 0.3, respectively. 

Fig. 4. Terminal velocity of the probe molecule as a function of the adsorption probability /. Nota- 
tions and values of g are the same as in Fig.3 except that the probe hopping probabilities are p = 0.98 
and q = 0.02. 

Fig. 5. Self-diffusion coefficient of the probe molecule as a function of the adsorption probability /. 
Notations and values of g are the same as in Figs. 3 and 4. 
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Fig.l. JSP, Benichou et al. 




FIG 2 JSP Benichou et al. 




FIG 3 JSP Benichou et al. 




FIG 4 JSP Benichou et al. 




FIG 5 JSP Benichou et al. 



